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Abstract. In this paper, the self-similar solution of resistive ad- 
vection dominated accretion flows (ADAF) in the presence of a pure 
azimuthal magnetic field is investigated. The mechanism of energy 
dissipation is assumed to be the viscosity and the magnetic diffusivity 
due to turbulence in the accretion flow. It is assumed that the mag- 
netic diffusivity and the kinematic viscosity are not constant and vary 
by position and a-prescription is used for them. In order to solve the 
integrated equations that govern the behavior of the accretion flow, a 
self-similar method is used. The solutions show that the structure of 
accretion flow depends on the magnetic field and the magnetic diffu- 
sivity. As, the radial infall velocity and the temperature of the flow 
increase, and the rotational velocity decreases. Also, the rotational ve- 
locity for all selected values of magnetic diffusivity and magnetic field 
is sub-Keplerian. The solutions show that there is a certain amount of 
magnetic field that the rotational velocity of the flow becomes zero. 
This amount of the magnetic field depends on the gas properties of the 
disc, such as adiabatic index and viscosity, magnetic diffusivity, and 
advection parameters. The solutions show the mass accretion rate in- 
creases by adding the magnetic diffusivity and in high magnetic pres- 
sure case, the ratio of the mass accretion rate to the Bondi accretion 
rate decreases as magnetic field increases. Also, the study of Lundquist 
and magnetic Reynolds numbers based on resistivity indicates that the 
linear growth of magnetorotational instability (MRI) of the flow de- 
creases by resistivity. This property is qualitatively consistent with 
resistive magnetohydrodynamics (MHD) simulations. 
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1. Introduction 

The standard geometrically thin, optically thick accretion disc model can success- 
fully explain most of the observational features in active galactic nuclei (AGN) and 
X-ray binaries (Shakura & Sunyaev 1973). In the standard thin model, the motion 
of the matter in the accretion disc is nearly Keplerian, and the gravitational energy 
released in the disc is radiated away locally. During recent years, another type of 
accretion flow has been studied in which it is assumed that energy released through 
dissipation processes in the disc may be trapped within accreting gas, and only the 
small fraction of the energy released in the accretion flow is radiated away due to in- 
efficient cooling, and most of the energy is stored in the accretion flow and advected 
to the central object. This kind of accretion flow is known as advection-dominated 
accretion flow (ADAF). The basis ideas of ADAF models have been developed by 
a number of researchers (e.g. Ichimaru 1977; Rees et al. 1982; Narayan & Yi 1994, 
1995; Abramowicz et al. 1995; Ogilvie 1999; Blandford & Begelman 1999). 

It seems that accretion discs, whether in star-forming regions, in X-ray bina- 
ries, in cataclysmic variables, or in the centers of active galactic nuclei, are likely 
to be threaded by magnetic fields. Consequently, the role of the magnetic fields on 
ADAFs has been analyzed in detail by a number of investigators (Bisnovatyi-kogan 
& Lovelace 2001; Akizuki & Fukue 2006, hereafter AF06; Shadmehri 2004, here- 
after Sh04; Ghanbari et al., 2007; Bu, Yung, & Xie 2009; Khesali & Faghei 2009, 
hereafter KF09). The existence of the toroidal magnetic field has been proven in the 
outer regions of the discs of young stellar objects (YSOs; Aitken et al. 1993; Wright 
et al. 1993) and in the Galactic center (Novak et al. 2003; Chuss et al. 2003). Thus, 
considering the accretion discs with a toroidal magnetic field have been studied by 
several authors (AF06; Begelman & Pringle 2007; Khesali & Faghei 2008, here- 
after KF08; Bu, Yung, & Xie 2009; KF09). 

The resistive diffusion of magnetic field is important in some systems, such as 
the protostellar discs (Stone et al. 2000; Fleming & Stone 2003), discs in dwarf 
nova systems (Gammie & Menou 1998), the discs around black holes (Kudoh & 
Kaburaki 1996), and Galactic center (Melia & Kowalenkov 2001; Kaburaki et al. 
2010). Also, two and three dimensional simulations of local shearing box have 
shown that resistive dissipation is one of the crucial processes that determines the 
saturation amplitude of the magnetorotational instability (MRI). As, linear growth 
rate of MRI can be reduced significantly because of the suppression by ohmic dis- 
sipation (Fleming et al. 2000; Masada & Sano 2008). 

AF06 proposed a self-similar advection-dominated accretion flow that the disc 
plasma is highly ionized, so they assumed that resistivity of the plasma is zero, 
and only viscosity is due to turbulence and dissipation in the disc. However, re- 
cent works represent importance of magnetic diffusivity in accretion discs (e.g. 
Kuwabara et al. 2000; Kaburaki 2000 ;Kaburaki et al. 2002; Sh04; Ghanbari 
et al. 2007; Krasnopolsky et al. 2010; Kaburaki et al. 2010). Sh04 studied a 
quasi-spherical accretion flow that dominant mechanism of energy dissipation was 
assumed to be the magnetic diffusivity due to turbulence in the accretion flow and 
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the viscosity of the fluid was completely neglected. The main focus of Sh04 was 
nonrotating accretion flow and ignored from toroidal magnetic field. Also, Sh04 
studied induction equation of magnetic field in a steady state that is not according 
to anti-dynamo theorem (e.g. Cowling 1981) and is usefull only in particular sys- 
tems where the magnetic dissipation time is very long, much longer than the age 
of the system. Ghanbari, Salehi, & Abbassi (2007) considered an axisymmetric, 
rotating, isothermal steady accretion flow, which contains a poloidal magnetic field 
of the central star and from toroidal magnetic field of the flow neglected. They as- 
sumed that the mechanisms of energy dissipation are the turbulence viscosity and 
magnetic diffusivity due to the magnetic field of the central star. They explored 
the effect of viscosity on a rotating disc in the presence of constant magnetic dif- 
fusivity. Similar to Sh04 they considered the flow in balance between escape and 
creation of the magnetic field, and ignored from toroidal component of magnetic 
field. 

As mentioned the observational evidences and the MHD simulations express 
that the toroidal component of magnetic field and the magnetic diffusivity are im- 
portant in accretion discs. Thus in this paper by using AF06 technique we will 
study the influence of presence of toroidal component of magnetic field in a vis- 
cous and resistive accreting gas, and investigate the role of non-constant magnetic 
diffusivity in systems that escaping and creating of magnetic field are not balanced. 
We will show that the present model from some aspects is qualitatively consistent 
with the observational evidences and the resistive MHD simulation results. The 
paper is organized as follows. In section 2, the basic equations of constructing a 
model for quasi-spherical magnetized advection dominated accretion flow will be 
defined. In section 3, self-similar method for solving equations which govern the 
behavior of the accreting gas was utilized. The summary of the model will appear 
in section 4. 

2. Basic Equations 

We use spherical coordinate (r, 9, tp) centered on the accreting object and make the 
following standard assumptions: 

(i) The flow is assumed to be steady and axisymmetric dt = = 0, so all flow 
variables are a function of only r ; 

(ii) The magnetic field has only an azimuthal component; 

(iii) The gravitational force on a fluid element is characterized by the Newtonian 
potential of a point mass, ^ = —GM*/r, with G representing the gravita- 
tional constant and M* standing for the mass of the central star; 

(iv) The equations written in spherical coordinates are considered in the equato- 
rial plane 9 = it/ 2 and terms with any 9 and 92 dependence are neglected 
(Ogilvie 1999; KF09). 
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(v) For the sake of simplicity, the self-gravity and general relativistic effects have 
been neglected; 

The behavior of such system can be analyzed by magnetohydrodynamics (MHD) 
equations. The general MHD equations are written as follows: 



! + V.(pv) = 0, 



<9v 
~dt 



+ (v • V)v 



1 



dB 

~dt 



V • B = 0, 



-Vp gas - pW + — J X B + F vis , 
J = V x B 
Vx(vxB-i)J), 

= Qdiss Qcool = fQ diss ^ 



(1) 

(2) 
(3) 
(4) 

(5) 
(6) 



where p, v, and p gas are the density, the velocity field, and the gas pressure, respec- 
tively; F v i s is the viscous force per unit volume; B is the magnetic field; J is the 
current density and 77 represents the magnetic diffusivity; 7 is the adiabatic index; 
The term on the right hand side of the energy equation, Qdiss, is the rate of heating 
of the gas by the dissipation and Q C ooi represents the energy loss through radiative 
cooling. The advection factor, / (0 < / < 1), describes the fraction of the dis- 
sipation energy which is stored in the accretion flow and advected into the central 
object rather than radiated away. The advection factor of / in general depends on 
the details of the heating and cooling mechanism and will vary with postion (e.g. 
Watari 2006, 2007; Sinha et al. 2009). However, we assume a constant / for sim- 
plicity. Clearly, the case / = 1 corresponds to the extreme limit of no radiative 
cooling and in the limit of efficient cooling, we have / = 0. 

Under the assumptions (i)-(v), the approximation of quasi-spherical symmetry, 
the equations (l)-(6) become 



r 2 dr 



(r 2 pv r 



dv r 



+ 



1 dp GM* 



dr p dr 



+ 



rQ 2 



— — — r (rB (p ), 
Airrp dr 



d 2 1 d 

pv r — {r il) = — — 

dr r z dr 



upr 



dr 



1 



7-1 

ld_ 

r dr 



dp 7p d 
r dr r 2 dr 



(rV) 



fQdissi 



rv r B v - r]—(rB v ) 



— B v . 



(7) 
(8) 
(9) 
(10) 
(11) 
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Here v r the radial velocity, p the mass-loss rate per unit volume, Q, the angu- 
lar velocity, the toroidal component of magnetic field, B v is the field escap- 
ing/creating rate due to a magnetic instability or dynamo effect, v the kinematic 
viscosity coefficient. We assume both of the kinematic coefficient of viscosity and 
the magnetic diffusivity due to turbulence in the accretion flow, it is reasonable to 
use these parameters in analogy to the a-prescription of Shakura & Sunyaev (1973) 
for the turbulent, 

v = a ^L(l + l3) 1 -^ (12) 

V = Vo^(l + ^ (13) 

where VLk = (GM^/r 3 ) 1 / 2 is the Keplerian angular velocity. Narayan & Yi (1995) 
used a similar form for kinematic coefficient of viscosity, i.e. v = a(p gas / p^lx), 
also Sh04 applied a similar form for magnetic diffusivity, i.e. rj = rio(p gas / pQx)- 
Thus in comparison to Narayan & Yi (1995) and Sh04 prescriptions, we are using 
the above equations for the kinematic coefficient of viscosity v and the magnetic 
diffusivity rj. The parameters of a and r? are assumed to be positive constants, and 
we assume that a, r]o < 1 (Campbell 1999; Kuwabara et al. 2000; Sh04; King et 
al. 2007). The parameter of f5[= p ma g/Pgas] is the degree of magnetic pressure, 
Pmag = B^/8it, to the gas pressure. Since, we will apply steady self-similar 
method for solving system equations, this parameter will be constant throughout the 
disc, but really this parameter varies by position (see, KF08, KF09). The parameter 
of p is a constant and states importance degree of total pressure in the kinematic 
viscosity and the magnetic diffusivity. Clearly, the case of p = 1 corresponds to 
traditional a-model and in case of p = total pressure is used. Note the kinematic 
coefficient of viscosity and the magnetic diffusivity are not costant and depend on 
the physical quantities of the flow, we will show the quantities of v and 77 in our 
self-similar solution vary with radius r 1 / 2 . 

The ratio of the kinematic coefficient of viscosity to the magnetic diffusivity 
is defined by the magnetic Prandtl number, P m = u/rj. By using equations (12) 
and (13) the magnetic Prandtl number in present model is P m = a/770 • We will 
consider conditions that P m = 00 (in case that magnetic diffusivity is zero), P m > 
1, and P m < 1. 

For the heating term in equation (10), Qdiss> we use two sources of dissipation: 
the viscous and resistive dissipations. Thus, for Qdiss we can write 

QMss = vpr 2 (^)\£.J 2 (14) 

where first term is related to viscous dissipation and second term is related to resis- 
tive dissipation. By converting the gas pressure and the magnetic pressure in terms 
of sound speed (c 2 = p ga s/p) and Alfven speed (c\ = B^/Anp), and by using 
equations (12)-(14) for the equations (7)-(10), we can write 
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dr J 
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Qdiss 



(1 + 0) 
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ar pc; 
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3. Self-Similar Solutions 

We seek self-similar solutions in the following forms: 



-c±a\ 





V r 




jGM* 


= c 3 


GM* 
r 


. B l _ 
47T/3 


2/3c 3 - 



(15) 
(16) 
(17) 
(18) 
(19) 

(20) 



(21) 

(22) 
(23) 
(24) 



where coefficients c±, C2, and C3 are determined later, we assume a power-law 
relation for density 



p(r) = Po r s 



(25) 



where po an d s are constant. By using above self-similar quantities, the mass-loss 
rate and the field escaping/creating rate must have the following form: 



P = Por s " 3 / 2 



(26) 
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B v = SorV (27) 

where po and Bq are constant. 

Substituting the above solutions in the continuity, momentum, angular momen- 
tum, energy, and induction equations [(15)-(20)], we can obtain the following rela- 
tions: 



Po = - + 2) apoCx^/GM*, 



(28) 



i c 2 a 2 + 1 _g + Ca ( 8 _ 1 + p( 1 + 8 ^ = 0) (29) 
ci = 3( S + 2)(l + /3) 1 -^c 3 , (30) 



-^ya Cl (2 S - 2 + 3 7 ) = 1/(1 + p) 1 ^ (dac 2 . + 2 % /3c 3 (l + s) 2 ) , (31) 

B = - S -GM* V2^pb~ 3 (2 Cl a + r/oc 3 (l + s)(l + Z?) 1 ^) . (32) 

Above equations express for s = —3/2, there is no mass loss, while for 
s > —3/2 mass loss (wind) exists. The escape and creation of magnetic fields 
are balanced in s = — 2 + (3??o)/(6a + 770) that Bq becomes zero, solving it in 
s = —3/2 (no wind) implies that 7/0 = (6/5)a. This amount of 770 corresponds 
to the magnetic Prandtl number 5/6. Thus, when the flow has s = —3/2 and 
Vo = (6/5)a, we expect the escape and creation of magnetic field are balanced and 
there is no mass loss. In 770 = for balance of escape and creation of magnetic 
field, we have s = —2. In this paper, only case of no wind (s = —3/2) is consid- 
ered that po = and B^ oc r -11 / 4 , we note B v in P rn = 5/6 will be zero, too. 
The equations (29)-(31) imply that for given a, 770, P, s, and / form a closed set of 
equations of c\, C2, and c 3 which will determine behavior of the accretion flow. 

3.1 Analysis 

By using the equations (29)-(31), the coefficients of Cj have the following forms 

ci = ^2 (d a + Jd 2 + 8a 2 ) , (33) 

c 2 = -^c 1 D 1 D 2 , (34) 

C3 = 57 1 ^\ ClDl ' (35) 
3(s + 2) 
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where 



A = (i + fl^ (36) 

(2g-2 + 3 7 ) r/o/3 (1 + s) 2 

° 2= (7-D/ (37) 

_ a -l + /3(l + a ) 

A3 ^ , (38) 

Da = %Vi (\d 2 + D 3 ) . (39) 



3 V3 

The obtained results imply that the model is parametrized by the ratio of specific 
heat, 7, the standard viscous parameter, a, the magnetic diffusivity parameter, r/o, 
the degree of magnetic pressure to gas pressure, (3, the advection parameter, /, and 
the mass-loss rate parameter, s. 

The equation (37) implies that the value of D% for a value of (3 will be zero. 
From equations (22) and (34) we can write ft oc c\ oc D 2 . Thus we conclude for a 
value of (3 that we define it fa (braking (3), the angular velocity will be zero. Over 
of /3b, c 2 , becomes negative that it is not physical. By solving D 2 = for the (3 
parameter, we can write 

3a (s + 2) ( 2s -2 + 3 7 ) 

A = "7^(- + i) 2 (7-i) ■ (40) 

In the case of no mass-loss, s = — 3/2, we can write 

& = ^fV^iV (41) 



/% V 7-1 

The above equation in terms of the magnetic Prandtl number becomes 

& = 18(%(^A (42) 



/ ' V 7-1 

The equations (40)-(42) express that the amount of the (3b parameter depends on 
/, P m , and 7. As /3b decreases by adding the advection degree, while increases by 
adding the magnetic Prandtl number. In a flow with high conductivity, that P m is 
very large, fib becomes very large. Also in the flow with high resistivity and low 
viscosity, f3b will be small. Since the magnetic pressure fraction always is positive 
or equal to zero ((3 > 0), so in presence model 7 < 5/3. 

Examples of the coefficients Cj in two cases of \i = and 1.0 are shown in 
figures (1) and (2) as a function of the degree of magnetic pressure to the gas 
pressure, (/?), for different value of the magnetic diffusivity, 770- In figures (1) and 
(2), C4 is the viscous torque that is obtained from right hand side of equation (30), 
C5 is total energy dissipation by the viscosity and the resistivity, and is calculated 
by right hand side of equation (31), and cq is the ratio of the resistive dissipation to 
the viscous dissipation. 
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3.1.1 First Case: /i = 

The parameter of \i has appeared in the equations of the kinematic coefficient of 
viscosity and the magnetic diffusivity to indicate important degree of total pres- 
sure in them. In the case of fi = 0, the kinematic coefficient of viscosity and the 
magnetic diffusivity become 

u = a-§-{l + ft) 
= ac 3 (l + /3) v / GM* r 1 / 2 (43) 

V = Vo^(l + ft) 
\Ik 

= r / oc 3 (l + /3)v / GM*r 1 / 2 . (44) 

In this case the function forms of the kinematic coefficient of viscosity and 
the magnetic diffusivity deviate by factor (1 + ft) from the function forms of used 
by NY95 and Sh04. Also the profiles of non-resistive and non-magnetic case is 
shown to compare the present model with canonical ADAF solutions (e.g. NY95). 
The existence of (1 + ft) in v and rj causes the viscosity and resistivity increase 
by adding the toroidal magnetic field. The v and r\ have direct effects on the vis- 
cous torque (C4) and the energy dissipation (C5). Thus, we expect increase of C4 
and C5 by adding the toroidal magnetic field that the profiles of C4 and C5 confirm 
it. Also, by adding the magnetic diffusivity (i]o) of the flow, C4 and C5 increase. 
Because the ohmic dissipation increases due to resistivity that it also increases the 
flow temperature (the profiles of C3 confirm it), and since the viscous torque is pro- 
portional with temperature (sound speed), thus the viscous torque, C4, increases by 
adding the magnetic diffusivity. The profiles of C4 and C5 imply that for all value 
of the ft and 770, the total energy dissipation and the viscous torque are larger than 
the canonical ADAF solutions. The profiles of C3 show that the temperature of the 
flow by adding the toroidal component of magnetic field decreases. This property 
is qualitatively consistent with the results of Bu et al. (2009) and KF09. The pro- 
files of cq show that the ratio of the resistive dissipation to the viscous dissipation 
increases by adding the toroidal component of magnetic field (ft) and the magnetic 
diffusivity (770). As in small amounts of magnetic field and the magnetic diffusiv- 
ity, the dominant heat generated is the viscous dissipation, while in large values of 
magnetic field and the magnetic diffusivity, the dominant heat generated will be the 
resistive dissipation. Also, figure (1) shows by adding the ft and t]q parameters, the 
radial infall velocity, c\, increases, and the angular velocity, C2, decreases. It can 
be due to increase of the viscous torque (C4) by parameters of ft and 770- The raise 
of the viscous torque by adding r/o and ft parameters generates a larger negative 
torque in angular momentum equation and causes the angular velocity of the flow 
decreases, and the matter accretes with larger speed. In the present model, the flow 
rotates slower than canonical ADAF and accretes speeder than it. The increase of 
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Figure 1. Physical quantities of the flow as a function of the degree of magnetic pressure 
to the gas pressure, for several values of ijo = 0, 0.12, 0.15, and 0.2 that corresponding to 
P m — oo, 5/6, 2/3, and 5/10. The disc density profile is set to be s = —3/2 (no wind), 
the ratio of the specific heats is set to be 7 = 1.3, the viscous parameter is a = 0.1, and 
the advection parameter is / = 1.0. 
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the radial infall velocity by adding the parameter of (3 is qualitatively consistent 
with the results by Bu et al. (2009) and KF09. 

3.1.2 Second Case: fx = 1 

In the case of \i = 1, for the kinematic coefficient of viscosity and the magnetic 
diffusvity we can write 



In this case the function forms of the kinematic coefficient of viscosity and the 
magnetic diffusivity are the same as Sh04 and NY95 used. The behavior of the 
physical quantities of the flow in this case are shown in figure (2). The absence of 
(1 + f3) in this case for v and rj in comparison to previous case causes the value 
of them decrease by factor (1 + f3). The effects of absence of this factor increases 
by adding parameter of (3. The profiles of the angular momentum transport (04) 
and total energy dissipation (05) imply that the amounts of them by factor (1 + (3) 
compared with previous case decrease. As their behavior in terms of the toroidal 
component of magnetic field have changed and they decrease by adding the (3 pa- 
rameter. However, the magnetic diffusivity has the previous effects and increase 
these two quantities. Also, the solutions show that the viscous torque and the to- 
tal dissipation in this case are smaller than canonical ADAF. The weakening of 
C4 and C5 by adding the (3 parameter reduces the radial infall velocity. Here, the 
radial infall velocity increases by adding 770 that is due to increase of C4 and C5. 
The behavior of other physical quantities in terms of the [3 and 770 parameters does 
not change, and represent small variations. To compare the radial infall velocity 
profiles with canonical ADAF solutions implies that the flow accretes slower than 
canonical ADAF that is different with previous case. The decrease of the tempera- 
ture and the viscous torque by adding the parameter of (3 is qualitatively consistent 
with the results of Bu et al. (2009). 

3.2 Mass Accretion Rate 

In according to assumptions of section (2) the mass accretion rate defines as 




ac 3y /GM* r 1 ' 2 . 



(45) 




(46) 



M = -4irr 2 pv r . 



(47) 



The mass accretion rate under self-similar transformations of equations (21) and 
(25) becomes 

M = Aitapocx y/GMl r s+3/2 . (48) 
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Figure 2. Physical quantities of the flow as a function of the degree of magnetic pressure 
to the gas pressure, for several values of ijo = 0, 0.12, 0.15, and 0.2 that corresponding to 
P m — oo, 5/6, 2/3, and 5/10. The disc density profile is set to be s = —3/2 (no wind), 
the ratio of the specific heats is set to be 7 = 1.3, the viscous parameter is a = 0.1, and 
the advection parameter is / = 1.0. 
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In our interesting case, s = —3/2 (no wind), for the mass accretion rate we can 
write 

M = A-TtapQCi^jGM* (49) 

that implies the mass accretion rate does not vary by position. This result is quali- 
tatively consistent with the results by Sh04, Ghanbari et al (2007), and AF06. Al- 
though the present model of accretion flow is different from Bondi (1952) accretion 
in various aspect, we can define Bondi accretion rate as 



M Bondi = ^Ml I P{Q0) 



cf(oo) 



5-37 



(5-3 7 )/2( 7 -l) 



(50) 



where p(oo) and c s (oo) are the density and the sound speed in the gas far away 
from the star (Frank et al. 2002). Bondi accretion rate in terms of our self-similar 
transformations becomes 



M Bondi = nVGM* 



po 

3/2 



-, (5-3 7 )/2( 7 -l) 



5-3 7 



(51) 



Thus, we can write the mass accretion rate in term of Bondi accretion rate as follows 



M/M B ondi = faciei 12 



5-37 



(3 7 -5)/2( 7 -l) 



,s+3/2 



In our interesting case, s = —3/2 (no wind), we can write 



c 7 = M/M Bo ndi = 4acic 3 7 



5 - 3 7 



(3 7 -5)/2( 7 -l) 



(52) 



(53) 



Here, we defined new parameter of c-j that indicates the ratio of the mass accretion 
rate to Bondi accretion rate. Examples of the coefficient of cj in two cases of p = 
and 1.0 are shown in figures (3) as a function of the degree of magnetic pressure to 
the gas pressure, (/3), for different value of the magnetic diffusivity, 770 . The profiles 
of cj show that in the case p = 0, the mass accretion rate increases by adding the 
toroidal magnetic field and the magnetic diffusivity. While the solutions for p = 1 
imply that the mass accretion rate decreases by adding the toroidal magnetic field 
and increases by adding the magnetic diffusivity. On the other hand, the magnetic 
diffusivity in two cases causes the mass accretion rate increase. In high magnetic 
pressure, the cj profiles for both cases show that the ratio of the mass accretion rate 
to the Bondi accretion rate is decreased with an increase in magnetic pressure. This 
property is qualitatively consistent with results of Kaburaki (2007). Comparision of 
the present model with canonical ADAF solutions implies that in the case of p = 
the flow accretes quicker than canonical ADAF, however in the case of p = 1 the 
flow accretes slower than canonical ADAF. Also, the profiles of c-j in both of cases 
show the mass accretion rate in our model is smaller than Bondi accretion rate that 
is in adapting with observational evidences from Sgr A*, M87, and NGC 4261 
(Kaburaki 2007). 
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Figure 3. The ratio of mass accretion rate to Bondi accretion rate (c-j = M j Msondd 
as a function of the degree of magnetic pressure to the gas pressure, for several values of 
770 = 0, 0.12, 0.15, and 0.2 that corresponding to P m = 00, 5/6, 2/3, and 5/10. The disc 
density profile is set to be s = —3/2 (no wind), the ratio of the specific heats is set to be 
7 = 1.3, the viscous parameter is a = 0.1, and the advection parameter is / = 1.0. 



In section 3.1, the upper limit of the magnetic field obtained and mentioned 
with fa. By substituting fa and equations (33)-(39) in equation (53) and assume of 
s = —3/2, the mass accretion rate to the Bondi accretion rate (07) approximately 
is 



<-■ M -V//;,,,,,/, ~ 21\'l2o <j\ ( ' ''" ''' 

= 24>/2 a 51 



vl-jU 



(5 + /3 6 ) 5/2 
(l+ 18 ^ p ") 1 " ft 



(5 + 18^1 



5/2 



where 



9i 



(3 7 -5) 

2 1 WTT) 



5-37 



92 



5/3-7 



7 



To obtain the root of derivative of c-j in terms of /, we can write 

18\ /3 + 2/i N 



df 



92 Pm, 



(54) 



5 J \l-2n, 

The above equation states the mass accretion rate in f max becomes maximum and 



Viscous and Resistive ADAF 



15 



for / > fmax the mass accretion rate decreases by advection degree. The f ma x 
for n > 1/2 becomes negative, while < / < 1. Thus, the relation of f max is 
valid only for the case of fi = 0. In the case of fi = 0, f max = (54/5) giP m 
and fa = 5/3. Thus, we expect in dominant magnetic case (/3 > 1), the accretion 
efficiency decreases by advection degree parameter. In the case of \i = 1, due to the 
lack of any extremum, the mass accretion rate only increases by advection degree 
and does not show any decrease. 

In the case of high magnetic pressure (/3 » 1), equation (54) becomes 

c 7 = M/M Bondi « 24\/2 a gi 

The above equation implies that the mass accretion rate to Bondi accretion is 
strongly depends on viscosity parameter, Prandtl number and the advection de- 
gree. The mass accretion rate increases by / and decreases by P rn . Thus, accretion 
efficiency increases by the advection degree parameter in high magnetic pressure. 




3.3 Timescales 

To estimate the effect of viscosity and resistivity on the accretion discs, we com- 
pare the viscous and resistive timescales with accretion timescale. The accretion 
timescale, t acc , and the viscous timescale, t„j SC , are given by 

r 

(55) 



— V r 



r 2 



tvisc — • (56) 
V 

We are using a similar functional form of t visc for the resistive timescale, t resis , 
that is given by 

tresis — • (57) 

V 

By using self-similar forms of physical quantities, we can write 

-(1 + /?) A 



tresis _^_J_^_, | , JV'-'-l 
tacc V0 C3 



= 3(-)(s + 2). (58) 
Vo 

The equation (30) is used for fraction of 01/03. As, we said in previous section, 
in present model a/770 is the magnetic Prandtl number, P rn , so above equation 
becomes 

t f^ = W m {s + 2). (59) 
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For our interesting case, s = —3/2 (no wind), we can write 



tresis 3 



-P m . (60) 



The above equation implies that for P m < 2/3, the magnetic diffusivity timescale 
is shorter than or equal to accretion timescale, while for P m > 2/3 the accretion 
timescale is shorter. 

Similar calculations for the viscous timescale express 

^=3( S + 2), (61) 

tacc 

where in no wind case (s = —3/2) becomes 

tirisf 



= (3/2). (62) 

L acc 

Thus, the viscosity timescale will be longer than the accretion timescale. To com- 
pare the magnetic diffusivity with the viscous timescales, we can write 

tresis T /fj 

tvisc T^/u 
V 

V 

= Pm- (63) 

Thus, the magnetic Prandtl number specifies which one is shorter. For exam- 
ple in flow with high conductivity (e.g. AF06; KF09), 77 — > 0, the magnetic 
Prandtl number limits to infinity, and so the magnetic diffusivity timescale will 
be very longer than the viscous timescale. On the other hand, for a flow with 
finite resistivity and tiny viscosity (e.g. Sh04), the magnetic Prandtl number lim- 
its to zero, and so the magnetic diffusivity timescale is very shorter than viscous 
timescale. When the resistivity and the viscosity are approximately equal, P m ~ 1, 
we expect 

tresis ^ t v i sc . Also, in special case of Pm — 5/6 and s — 3/2 
that escape and creation of magnetic field are balanced and there is no mass-loss, 

tresis = (5/6)t^j sc . 



4. Summary and Discussion 

In this paper, the influences of the resistivity on the structure of the advection- 
dominated accretion flow is investigated. It is used only azimuthal component 
of magnetic field that is consistent with observational evidence of Galactic center 
(Novak et al 2003; Chuss et al. 2003; Yuan 2006). The a-prescription is used for 
the kinematic coefficient of viscosity and the magnetic diffusivity. The equations of 
the model are solved by a semi-analytical self-similar method in comparison with 
the self-similar solution by AF06. 
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The physical quantities of disc are sensitive to the amounts of the magnetic 
pressure fraction (/?) and the magnetic diffusivity (r/ ) parameters. As, the angular 
velocity of the flow by adding the (3 and t]q parameters decreases. For a value of the 
magnetic pressure fraction, the angular velocity of disc becomes zero. This amount 
of the magnetic pressure fraction strongly depends on the properties of the accret- 
ing gas, such as the viscosity, resistivity, adiabatic index, and advection degree. 
The solutions represent the radial infall velocity increases by adding the magnetic 
diffusivity. Also the solutions show that the temperature of the flow decrease by 
adding the toroidal component of magnetic field. This result qualitatively is consis- 
tent with the results of Bu et al. (2009) and KF09. The profiles of the temperature 
of the flow show that it increases by adding the magnetic diffusivity that is due to 
the raise of the resistive dissipation. Comparison of the present model with Bondi 
accretion implies that for all values of the j3 and r] parameters, the Bondi accre- 
tion rate is larger than the mass accretion rate that is in accord with observational 
evidences of Sgr A*, M87, and NGC 4261 (Kaburaki 2007). Also, the mass ac- 
cretion rate profiles at high magnetic field express that the magnetic field reduces 
the mass accretion rate that is similar to results of Kaburaki (2007). We found that 
in the small magnetic field, the more heat generated in the flow is due to the vis- 
cous dissipation, while the ohmic dissipation will be dominant in large amounts of 
magnetic field and resistivity. 

As noted in the introduction, the MHD simulations show that linear growth of 
MRI decreases significantly by ohmic dissipation, linear growth of the MRI in the 
resistive fluid can be characterized by the Lundquist number (Smri = c\/r]Q) 
and magnetic Reynolds number {Reu = c 2 s /i]^l), where ca, c s , t], and Q have 
usual meaning. In terms of our self-similar transformations, the Lundquist num- 
ber and magnetic Reynolds number become Smri = 2/?/?7oC2(l + and 
Reu = 1/%C2(1 + M - The solutions of present model show that Smri and 
ReM decrease by resistivity. This property is qualitatively consistent with MHD 
simulation results (Fleming et al. 2000; Masada & Sano 2008). 

Here, latitudinal dependence of physical quantities is ignored, while some au- 
thors showed that latitudinal dependence is important in structure consideration of 
a disc (Narayan & Yi 1995; Sh04; Ghanbari et al. 2007). One can investigate 
latitudinal behavior of such discs. Furthermore, in a realistic model the advection 
parameter / is a function of position, one can consider such discs. 



Acknowledgements 

I would like to thank the referee for very useful comments that helped me to im- 
prove the initial version of the paper. I would also like to thank Markus Flaig for 
helpful discussion. 



18 



K. Faghei 



References 

Abramowicz, M., Chen, X., Kato, S., Lasota, J. P., Regev, O., 1995, ApJ, 438, L37 

Aitken D. K., Wright C. M., Smith C. H., Roche P. R, 1993, MNRAS, 262, 456 

Akizuki, C, Fukue, J., 2006, PASJ, 58, 469 (AF06) 

Balbus S. A., Hawley J. F., 1992, ApJ, 392, 662 

Begelman, M. C, Pringle, J.E., 2007, MNRAS, 375, 1070 

Bisnovatyi-Kogan,G. S., Lovelace, R. V. E., 2001, New Astron. Rev. 45, 663 

Blandford, R. D., Begelman, M. C. 1999, MNRAS, 303, LI 

Bondi, H., 1952. MNRAS, 112, 195 

Bu D., Yuan F, Xie F, 2009, MNRAS, 392, 325 

Campbell C. G., 1999, Geophys. Astrophys. Fluid Dynamics, 90, 113 

Chuss, D. T., Davidson, J. A., Dotson, J. L., Dowell, C. D., Hildebrand, R. H., Novak, G, 

Vaillancourt, J. E., 2003, ApJ, 599, 1 1 16 
Cowling, T. G, 1981, ARA&A, 19, 115 
Fleming T. P., Stone J. M., Hawley J. F, 2000, ApJ, 530, 464 
Fleming T. P., Stone J. M., 2003, ApJ, 585, 908 

Frank J., King A. R., Raine D., 2002, Accretion Power in Astrophysics. Cambridge Univ. 

Press, Cambridge 
Gammie C. F, Menou K., 1998, ApJ, 332, 659 
Ghanbari J., Salehi F, Abbassi S., 2007, MNRAS, 381, 159 
Ichimaru, S. 1977, ApJ, 214, 840 
Kaburaki, O., 2000, ApJ, 531, 210 
Kaburaki O., 2007, ApJ, 662, 102 

Kaburaki, O., Yamazaki,N., Okuyama, Y. 2002, NewA, 7, 283 
Kaburaki O., Nankou T., Tamura N., Wajima K., 2010, PASJ, 62, 1 177 
Khesali, A., Faghei, K, 2008, MNRAS, 389, 1218 (KF08) 
Khesali, A., Faghei, K., 2009, MNRAS, 398, 1361 (KF09) 
King A. R., Pringle J. E., Livio M., 2007, MNRAS, 376, 1790 
Krasnopolsky R., Li Z. Y, Shang H., 2010, ApJ, 716, 1541 
Kudoh T., Kaburaki O., 1996, ApJ, 460, 199 

Kuwabara T., Shibata K, Kudoh T., Matsumoto R., 2000, PASJ, 52, 1 109 
Masada Y, Sano T, 2008, ApJ, 689, 1234 
Narayan, R., Yi, I., 1994, ApJ, 428, L13 
Narayan, R., Yi, I., 1995, ApJ, 453, 710 

Novak, G., Chuss, D. T., Renbarger, T., et al. 2003, ApJL, 583, L83 
Ogilvie, G. 1., 1999, MNRAS, 306, L90 

Rees, M. J., Begelman, M. C, Blandford, R. D., Phinney, E. S. 1982, Nature, 295, 17 

Shadmehri, M., 2004, A&A, 424, 379 (Sh04) 

Shakura, N.I., Sunyaev, R.A., 1973, A&A, 24, 337 

SinhaM., Rajesh S. R., Mukhopadhyay B., 2009, RAA, 9, 1331 

Stone J. M., Gammie C. F, Balbus S. B., Hawley J. F, 2000, in Protostars and Planets IV, 

ed, V. Manning, A. Boss, S. Russell, Tuncsun: Uni. Arizona Press 
Watarai, K. 2006, ApJ, 648, 523 
Watarai, K. Y. 2007, PASJ, 59, 443 

Wright C. M., Aitken D. K, Smith C. H., Roche P. F, 1993, PASA, 10, 247 
Yuan F, 2006, JPhCS, 54, 427 



